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Water is a life-giving source, fundamental to human existence, yet, over a billion people 
lack access to clean drinking water. Present techniques for water treatment such as piped, 
treated water rely on time and resource intensive centralized solutions. In this work, we 
propose a decentralized device concept that can utilize sunlight to split water into hydrogen 
and hydrogen peroxide. The hydrogen peroxide can oxidize organics while the hydrogen bub¬ 
bles out. In enabling this device, we require an electrocatalyst that can oxidize water while 
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suppressing the thermodynamically favored oxygen evolution and form hydrogen peroxide. 
Using density functional theory calculations, we show that the free energy of adsorbed OH* 
can be used as a descriptor to screen for selectivity trends between the 2 e~ water oxidation to 
H2O2 and the 4 e“ oxidation to O2. We show that materials that bind oxygen intermediates 
sufficiently weakly, such as Sn 02 , can activate hydrogen peroxide evolution. We present a 
rational design principle for the selectivity in electrochemical water oxidation and identify 
new material candidates that could perform H2O2 evolution selectively. 

The global energy consumption is projected to increase with the increased energy con¬ 
sumption being concentrated in areas that rank highest on the water scarcity index. A key 
challenge in providing clean drinking water is to hnd a low-cost, energy-efficient approach to 
treatment and disinfection of water so that it is suitable for consumption. Providing piped, 
treated water requires time and resource intensive centralized facilities and an infrastructure 
that does not exist in many places today. ^ Conventional techniques for water disinfection 
typically involve the use of chlorine or ozone as the oxidant.^ However, chlorination gener¬ 
ates hazardous and carcinogenic compounds,^ while the use of ozone, though efficient and 
harmless, is expensive.^ Hydrogen peroxide is an attractive candidate for water treatment 
as its degradation product is water and it has strong oxidative properties. ® 

H2O2 is generated at an industrial scale through the oxidation of anthraquinone. This 
process is a multi-step method requiring signihcant energy input and generates substantial 
waste and its transportation causes possible hazards. ® A direct efficient and economic route 
for production of hydrogen peroxide could enable its wide-spread application in water treat¬ 
ment and other areas. The direct thermal catalytic synthesis from H2 and O2 on palladium 
based materials have been studied for many years. However, selectivity and production 
rate of H2O2 are far below the desired limit.® An alternate synthetic route is through direct 
electrochemical reduction of oxygen and protons.®’^® Experimental and theoretical studies 
have shown it is possible to selectively activate hydrogen peroxide generation from oxygen 
reduction, however, this route requires electrocatalysts that are made out of expensive metals 
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Figure 1 : Schematic of a novel ‘dream’ device concept that can absorb photons from sunlight 
and use it to split water into hydrogen and hydrogen peroxide. Hydrogen peroxide will 
decompose the organics in water and thus cleaning water and hydrogen will bubble out. 

such as gold, platinum and palladium. 

Our core idea is to identify a synthetic route for hydrogen peroxide that uses water as its 
only reactant and nses sustainable electricity preferably from sunlight. In Fignre 1 , we show 
the schematic of a decentralized water treatment device. This device employs a material 
that can absorb photons and generate electrons and holes with appropriate energy snch that 
the electrons can reduce protons to hydrogen and holes can oxidize water to H2O2. To 
enable this device, two material challenges need to be overcome. The hrst is we require a 
photon absorber whose band positions are suitably aligned such that it can catalyze hydrogen 
evolution and hydrogen peroxide evolution. The second, more formidable requirement is of 
an electrocatalyst that can catalyze H2O2 evolution and suppress the thermodynamically 
favored O2 evolution. The second electrocatalyst requirement forms the focns of our present 
work. 

In this work, using thermodynamic analysis based on density functional theory calcula¬ 
tions, we demonstrate the existence of material candidates that can activate H2O2 evolution 
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through the oxidation of water. We show that the free energy of adsorbed OH* can be used 
as a descriptor, to a hrst approximation, for determining trends in 2e“ vs 4e“ oxidation of 
H 2 O. We identify materials that are good candidate materials for H 2 O 2 generation. Among 
these, we identify Sn 02 and Ti 02 as candidate materials with good selectivity and our 
analysis provides a quantitative foundation for the identihcation of more efficient, selective 
electrocatalyst materials. This analysis provides a necessary, but not sufficient criterion for a 
good selective electrocatalyst. Undoubtedly, kinetic barriers are important for determining 
selectivity, however, such a thermodynamic analysis has proved successful in determining 
selectivity trends for oxygen reduction. 


Results 

An ideal electrocatalyst for the 4 electron oxygen evolution reaction should be capable of 
facilitating oxidation of H 2 O just above the equilibrium potential of 1.23 V. As a mini¬ 
mum requirement, the four charge transfer steps should have reaction free energies of the 
same magnitude equal to the equilibrium potential of 1.23 eV. We consider the associative 
mechanism shown belowU®^^® 


2H20(1) + * ^ OH* + H20(1) + H+ + e”, (la) 

OH* + H20(1) + H+ + e" ^ O* + H20(1) + 2H+ + 26", (lb) 

O* + H20(1) + 2H+ + 2e- ^ OOH* + 3H+ + de”, (Ic) 

OOH* + 3H+ + 3e- ^ 02 (g) + 4H+ + de” + *. (Id) 


However, an ideal electrocatalyst for the two electron oxidation of water to hydrogen peroxide 
should facilitate the oxidation just above the equilibrium potential of 1.77 V. This implies 
that each of the two charge transfer steps must have a reaction free energy of 1.77 eV. We 
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consider a similar associative mechanism for H 2 O 2 production, 


2H20(1) + * ^ OH* + H20(1) + H+ + e”, 


(2a) 


OH* + H20(1) + H+ + e“ ^ H202(1) + 2H+ + 2e“. (2b) 
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Figure 2: Free energy diagram of water oxidation plotted at U = 1.77 V, versus the reversible 
hydrogen electrode on SnO2(110) and RuO2(110). On Sn02(lll), the limiting step for 2e“ 
oxidation is the activation of H 2 O as OH* while that for the 4e“ oxidation is the oxidation 
of OH* to O*. On Ru 02(111), the limiting step for 2e“ oxidation is the formation of H 2 O 2 
from OH* while that for the 4e“ oxidation is the oxidation of O* to OOH*. 


In Figure 2, we show the calculated free energy diagram for the 2e“ and 4e“ oxidation 
of water on rutile type Sn02 and Ru02 at U = 1.77 V, versus the reversible hydrogen 
electrode. Using the free energy diagrams, an important parameter, the thermodynamic 
limiting potential, Ul, can be extracted and this is dehned as the lowest potential at which 
all of the reaction steps are downhill in free energy. This approach has been successfully used 
to rationalize trends in hydrogen evolution, oxygen reduction and oxygen evolution. ^'^> 20-22 
The potential determining step for the 2e“ oxidation on Sn02 is the activation of H 2 O as 
OH* while that for the 4e“ oxidation is the oxidation of OH* to O*. It is to be noted that the 
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calculated limiting potential for the 2e“ oxidation is lower than that for the 4e“ oxidation 
and we would expect this to show selectivity towards hydrogen peroxide generation. In 
contrast, the potential determining step for 2e“ oxidation on Ru 02 is the formation of H 2 O 2 
from OH* while the potential determining step for 4e“ oxidation is the oxidation of O* to 
OOH*. For RuO 2(110), the activation to OH* is facile at 1.77 V, however, because the 
further oxidation of OH* to O* is strongly exothermic, selectivity for H 2 O 2 is expected to 
be low on Ru02. Materials with an OH* binding energy between that on Sn02 and RUO 2 
are expected to have improved activity for the 2e“ oxidation provided selectivity for the 4e“ 
oxidation can be suppressed. 

The trends for oxygen electrochemistry is determined by the binding of three key in¬ 
termediates, O*, OH* and OOH*.^® However, it has been demonstrated that the binding 
of these intermediates on oxide materials are correlated. This enables the activity, given 
by the limiting potential, Ul, to be plotted as a function of a single descriptor, to a hrst 
approximation. These plots allow for the quantitative determination of the descriptor values 
that yield optimal catalyst activity. 

Generalizing this analysis, for the 4e“ oxidation of water, in the case of materials that 
bind oxygen intermediates too strongly, we have step Ic associated with the oxidation of 
adsorbed O* being the limiting step. Therefore, the free energy difference of the limiting 
step is given by, 

AGic = AGooh* — AGo* • (3) 

In the case of the materials that bind oxygen intermediates too weakly, we have step lb 
associated with the oxidation of adsorbed OH* being the limiting step for the 4e“ oxidation 
of water. Therefore, the free energy difference of the limiting step is given by, 

AGib = AGo* — AGoh* (4) 

For the 2e“ oxidation of water, the activity of materials that bind oxygen intermediates 
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too strongly, reaction 2b associated with the oxidation of OH* to H 2 O 2 is the limiting step. 
The free energy of the limiting step is given by, 

^G2h = ^Oh202(i) ~ AOoh* • (5) 

The activity of materials on the weak binding leg of the volcano is limited by 2a associated 
with the activation of water as OH* and the free energy difference of the limiting step is 
given by, 

AG'2a = AGqh* ~ (6) 

In this work, we have chosen free energy of OH*, AGoh*, as the descriptor and this choice 
is made as the activity for hydrogen peroxide evolution is directly determined by AGqh*- 
In Figure 3, we show the calculated thermodynamic limiting potentials, Ul, as a function of 
the free energy of OH*, AGoh*- 

The calculated activity volcano shows that the optimal catalyst for peroxide evolution 
should exhibit a free energy of OH* adsorption of 1.77 eV. We identify RUO 2 , Pt02 as mate¬ 
rials that bind oxygen too strongly for peroxide evolution. Pt 02 calculations are performed 
on a rutile structure. According to our analysis, we also identify Sn02 and Ti02 as materials 
that bind oxygen too weakly and have slightly lower potentials for peroxide evolution than 
oxygen evolution. It is worth highlighting that kinetic barriers are also important for deter¬ 
mining selectivity and the activity volcano analysis presented here provides a necessary, but 
not sufficient criterion for a good selective electro catalyst. 

Discussion 

The only direct experimental demonstration for selectively making H 2 O 2 over O 2 through 
electrochemical oxidation of water has been the recent work of MacFarlane and co-workers 
who employed MnOa, electrode with an ionic liquid based electrolyte to tune the thermo- 


7 




Figure 3: Activity volcano for the 2e and 4e oxidation of water. The scaling relations 
nsed to construct the volcano are AGooh = AGqh + 3.2 [eV] and AGq = 2AGoh + 0.28 
[eV].24 

dynamics and showed H 2 O 2 formation. Based on onr analysis, we hnd that the limiting 
potentials for hydrogen peroxide evolution and oxygen evolution on Mn 02 are quite close 
and small changes in surface energetics due to solvent effects (ionic liquid) could affect the 
selectivity. There has been no direct experimental demonstration of electrochemical H 2 O 2 
generation in a purely aqueous system. It is worth pointing out that photogeneration of 
H 2 O 2 has been demonstrated on Ti 02 and ZnO.^^’^® 

It has also been demonstrated that electrodes that have high overpotentials for oxygen 
evolution exhibit enhanced performance towards the decomposition of organics although the 
exact details of the mechanism are still unclear. For instance, it has been shown that 
Sn 02 which is poor at catalyzing oxygen evolution exhibits a higher efficiency for organic 
degradation compared to Pt.^° We suggest that generation of hydrogen peroxide is a pre¬ 
cursor to the decomposition of organics and based on this assumption, our analysis attributes 



the enhanced efficiency for the decomposition of organics to the weaker binding of oxygen 
intermediates on Sn02 compared to Pt. As a result, we would expect little decomposition 
of organics on Pt 02 and signihcant decomposition on Sn 02 . 

Our descriptor based approach can be used to identify possible candidate materials that 
could be effective at selectively catalyzing peroxide evolution. An avenue could be the doping 
of metal ions on cheap material such as Ti 02 , which has been pursued for oxygen evolution. 
We have searched for suitable doped Ti02 candidates for selective peroxide evolution. Based 
on our analysis, we identify Ti02 doped with Ru or Ir are interesting candidate materials that 
could exhibit enhanced activity for peroxide evolution, as discussed in the SI. In addition to 
selectivity, stability is a stringent requirement for water oxidation electrocatalystsand we 
expect the stability requirement to be play a crucial role in the selection of an electrocatalyst. 

Finally, the provided device requires a photon absorber whose band positions are suitably 
aligned for hydrogen peroxide generation. This poses a requirement for a photon absorber 
that the valence band maximum is greater than 1.77 V vs normal hydrogen electode (NHE) 
and the conduction band minimum is less than 0 V vs NHE. There are many material 
candidates that satisfy this criterion, including Sn02. Hence, identifying a suitable photon 
absorber is less challenging than hnding a selective electrocatalyst. This analysis suggests 
that it is possible to identify a single material that can carry out the catalysis as well as the 
photon absorption. 

We have outlined a quantitative framework for determining selectivity in electrochemical 
water oxidation. We show that it is possible to selectively catalyze the 2e“ oxidation to 
H 2 O 2 over the thermodynamically favored oxygen evolution. This can be accomplished under 
certain range of potentials by choosing catalysts that are inefficient at carrying out oxygen 
evolution. We show Sn02 and Ti02 as materials that exhibit suitable bonding characteristics 
for peroxide evolution and identify doped Ti 02 candidate materials that could carry out this 
process more efficiently. This shows that it is possible to selectively form fuels or chemicals 
that involve a smaller number of proton-coupled electron transfer over its thermodynamically 
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favored competing reaction that involves a larger number of proton-electron transfer. We 
expect this core idea to be broadly useful given the ubiquity of adsorbate scaling relations 
and we expect it to be particularly useful for nitrogen and carbon electrochemistry. 


Methods 

Free energy diagrams: The free energy diagram is calculated based on density functional 
theory calculations which accounts for zero point energy and entropic corrections and a 
detailed discussion is presented in the SI. The effect of potential, U, is included by shifting 
the free energy of an electron by -eU and the free energy at a potential U, can thus be 
calculated using the relation, AG = AG° — eU where U is the potential relative to the 
reversible hydrogen electrode and is the calculated reaction free energy under standard 
conditions. 

Computational details: A detailed description of the computational details is given 
in the SI. 
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Computational Details 

Formation energies of OH* and O* on SnO2(110) and OH*, O* and OOH* MnO2(110) are 
obtained from Density Functional theory calculations by Man et al. ^ Additional calculations 
are performed on SnO2(110) to obtain the formation energy of OH* . It is found that OOH* 
spontaneously transfers a hydrogen atom to a nearby bridging oxygen atom. The free energy 
of OOH * on SnO2(110) is therefore approximated from the scaling relation between OH * 
and OOH * and the formation energy of OH * 

The computational details used are similar to those in ref. 1 and given below for com¬ 
pleteness. 

Density functional theory calculations are performed with the DACAPO DFT code. Ionic 
cores are described using Vanderbilt ultrasoft pseudo potentials and Kohn-Sham states are 
expanded in plane waves with an energy cutoff of 350 eV, while the electron density is 
expanded in plane waves with an energy cutoff corresponding to 500 eV.^ Occupation of 
one-electron states follows a Fermi-Dirac distribution with ksT = 0.1 eV and total energies 
are extrapolated to ksT = OeV. Effects of exchange and correlation are described using the 
RPBE functional."^ 

The SnO2(110) surface is modeled using a slab with a (1 x 2) surface supercell cell 
consisting of 4 trilayers. The geometry of the bottom two trilayers is hxed in the bulk 
position, and adsorbates are added to the topside of the slab. Slabs are separated by 15 A 
of vacuum and the electrostatic dipole interaction between periodically repeated slabs has 
been removed.^ The hrst Brilluin zone is sampled using 4x4x1 Monkhorst-Pack grid of 
k-points.® Adsorbates and the two topmost trilayers are optimized until the maximum force 
component is below 0.05 eV/A. 
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Free Energy Diagrams 

Potential dependent free energies are calculated using the computational hydrogen electrode 
reference. ^ Consider the initial step in water oxidation, written here in acid 

H20(1) + * -^ OH * + H+ + e”. (1) 

At 0 V vs an reversible hydrogen electrode (RHE) the reaction 

iHj(g)-> H+ + e- (2) 

is in equilibrium, so the chemical potential of a proton and an electron is equal to the 
chemical potential of 5 H 2 (g) 

= hH+ +he-- (3) 

Therefore the reaction free energy of eq. 1 at 0 V versus RHE, AG°, can be calculated from 
the equivalent gas phase reaction 

H20(1)+ * - ^OH* +|H2(g) (4) 

AG° = G(OH * ) — G(*) — MHaOO) + 2^H2(g)- (5) 

At an arbitrary potential, U versus RHE, the chemical potential of an electron is shifted by 
—ell and correspondingly the reaction free energy of eq. 1 is given by 

AG = AG^ - eU. (6) 

The approach is easily generalized to include formation of O* and OOH* adsorbates. 

The reaction free energy of eq. 4 is calculated from DFT simulations by adding con¬ 
tributions from entropy and vibrational zero-point energies (ZPE) to the reaction energies 
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obtained from DFT. 


— /S.Edpx Aii/zpE — T/\S. 


(7) 


For adsorbates, the zero point energy is calculated from vibrational frequencies calculated 

on RuO 2 ( 110 ) and taken from ref. 1, while the entropy of adsorbed species is assumed to 

be negligible. For molecules, the ZPE is obtained from DFT calculated frequencies, while 

the entropy is taken from experiment.® The chemical potential of H20(l) is calculated as 

the chemical potential of H20(g) at 0.035 bar, which is the vapor pressure of H 2 O at room 

temperature. The contributions from ZPE and entropy are listed in SI. 

Table SI: Zero point energies and entropic contributions to adsorbates and water 
oxidation reactions. Energies are in eV. 



TS 

T^S 

Ezpe 

^Ezpe 

AEzpe - TAS 

H20(1) 

0.67 


0.56 



H2(g) 

0.41 


0.27 



OH* 

0.0 


0.36 



0 * 

0.0 


0.07 



OOH* 

0.0 


0.43 



OH* +|H2(g) 

0.20 

-0.47 

0.50 

-0.06 

0.41 

0* +H2(g) 

0.41 

-0.27 

0.34 

-0.22 

0.05 

OOH* +iH2(g) 

0.61 

-0.73 

0.83 

-0.29 

0.45 
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Screening for new candidate materials 

Based on the activity volcano, we can screen for new candidate materials that could posses 
higher selectivity than Sn02. One strategy is to use a cheap material such as Ti02 and 
using doping to tune the adsorption energy. Our analysis suggests that a strengthening of 
the OH* binding energy by about ~0.3 eV relative to that of Ti02 could lead to substantially 
improved selectivity and electrocatalytic activity for H 2 O 2 evolution. 

We screened the calculated adsorption energies on doped rutile TiO2(110) surfaces based 
on a substitutional model with 6.25% transition-metal impurities relative to the host Ti 


4 



atoms in the slab, that is, M-Tii5032.® The analysis considered transition metals M=V, Nb, 
Ta, Cr, Mo, W, Mn, Fe, Ru, Ir, and Ni as dopants at four different substitutional sites. 
Among the considered cases, our analysis suggests the most promising candidates are Ir or 
Ru-doped Ti02 as shown in Figure SI. 



Figure SI: Activity volcano for the 2e“ and 4e“ oxidation of water with the identihed doped 
Ti02 candidates. The adsorption site for the doped Ti02 candidates, Ir and Ru, is on top 
of a 5-fold coordinated Ti site while the Fe-doped case is on top of a 6-fold coordinated Ti 
site. 
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